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ABSTRACT 

Classical nova events in symbiotic stars, although rare, offer a unique opportunity to probe the 
interaction between ejecta and a dense environment in stellar explosions. In this work, we use X-ray 
data obtained with Swift and Suzaku during the recent classical nova outburst in V407 Cyg to explore 
such an interaction. We find evidence of both equilibrium and non-equilibrium ionization plasmas 
at the time of peak X-ray brightness, indicating a strong asymmetry in the density of the emitting 
region. Comparing a simple model to the data, we find that the X-ray evolution is broadly consistent 
with nova ejecta driving a forward shock into the dense wind of the Mira companion. We detect a 
highly absorbed soft X-ray component in the spectrum during the first 50 days of the outburst that 
is consistent with supersoft emission from the nuclear burning white dwarf. The high temperature 
and short turn off time of this emission component, in addition to the observed breaks in the optical 
and UV lightcurves, indicate that the white dwarf in the binary is extremely massive. Finally, we 
explore the connections between the X-ray and GeV 7-ray evolution, and propose that the gamma 
ray turn-off is due to the stalling of the forward shock as the ejecta reach the red giant surface. 

Subject headings: stars: white dwarfs, X-rays: stars, ultraviolet: stars 


1. INTRODUCTION 

Classical novae (CNe) are a subclass of cataclysmic 
variable star (white dwarfs accreting material from a bi- 
nary companion) that experience rare, sudden increases 
in brightness of up to 9 magnitudes. These outbursts are 
now understood to be the result of a thermonuclear run- 
away event in the shell of accreted material on the white 
dwarf surface, initiated once the temperature and pres- 
sure required for nuclear burning are reached. A small 
number of nova systems have been detected in outburst 
more than once, and are referred to as recurrent novae. 
However, all novae are expected to recur on timescales 
ranging from a less than a decade to 100,000s of years, de- 
pending primarily on the mass of the white dwarf ( Mwd ) 
and the mass accretion rate (M). More massive white 
dwarfs have higher surface gravities, and therefore higher 
pressures in the accreted shell - as a result, less mass is 
needed to reach the critical pressure required for the on- 
set of the thermonuclear runaway. Higher mass accretion 
rates decrease the time it takes to reach the critical ac- 
creted mass. 

Most classical nova events occur in close binary sys- 
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terns, where the mass donor is a Roche Lobe filling, late 
type main sequence star. More rarely, the system is a 
wide binary hosting a white dwarf and a late type giant 
star, which may or may not fill its Roche Lobe. Four 
of the known recurrent novae (RS Oph, T CrB, V745 
Sco and V3890 Sgr) have giant companions. These sys- 
tems are also classified as symbiotic stars. The nature 
of the mass donor has major implications for the evolu- 
tion of the outburst over time, particularly in terms of 
the interaction of the nova ejecta and the circumbinary 
environment. 

Over the last three decades, CNe in outburst have been 
established as an important class of X-ray emitting ob- 
ject. They are perhaps best known as supersoft X-ray 
sources (SSS), characterized by blackbody like emission 
with 40 < kT (eV) < 80 and 10 36 < Lx (erg s _1 ) < 
10 38 . The SSS phase of the nova outburst occurs when 
the accreted envelope has expanded enough to become 
optically thin to soft X-rays, and thus we see the ra- 
diation originating in the burning layer. CNe are also 
known to emit hard X-rays at early times, typically at 
much lower luminosities than associated with the SSS 
phase (Mukai & Ishida 2001; Mukai et al. 2008). This 
emission has two possible origins. In most CNe, where 
there is very little material surrounding the binary, the 
X-ray emission is thought to originate in internal shocks 
within the nova ejecta. More rarely, the X-rays are 
due to interaction of the nova ejecta with a dense cir- 
cumbinary medium, as in RS Oph (Sokoloski et al. 2006; 
Bode et al. 2006). Events of this type can be considered 
analogous to supernovae in dense circumstellar media, 
whose early X-ray emission is also due to the interaction 
of the ejecta and the local environment (Chevalier 1982; 
Fransson et al. 1996). 

Novae in symbiotic star systems offer a unique opportu- 
nity to study shock physics relevant for supernova rem- 
nant evolution. The local nature of novae make them 
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easier to study in X-rays, and the much smaller ejected 
masses involved allow the evolving interaction with the 
surrounding medium to be observed over modest time pe- 
riods. Symbiotic stars may also be progenitors of some 
type la supernovae. In recent years, variable Na D ab- 
sorption features have been detected in a small number of 
systems that have been interpreted as the presence of ma- 
terial around the progenitor system, (see e.g. Patat et al. 
2007; Simon et al. 2009). Patat et al. (2011) presented 
similar evidence for evolving Na D absorption in the 2006 
outburst of RS Oph, and discussed the role of symbiotic 
novae as the progenitors of at least a subset of SNe la. 

1.1. The 2010 outburst of V407 Cyg 

V407 Cyg is a member of the small group of symbiotic 
Mira stars (binaries comprised of a white dwarf and a 
Mira- type AGB companion). The system is only weakly 
symbiotic, with no detectable continuum from the hot 
component in quiescence in the optical or near UV bands, 
and only weak Balmer series emission lines. The He II 
4686 emission line is occasionally detected. The pulsa- 
tion period of the Mira is 745 ± 15 days (Munari et al. 
1990; Kolotilov et al. 1998). This is the longest pulsation 
period of all known symbiotic Mira systems. Lithium 
has been identified in the optical spectrum of V407 Cyg 
(Tatarnikova et al. 2003a), indicating that this system 
may be experiencing so-called “hot bottom burning” 
(HBB.) This makes determining the distance to the sys- 
tem difficult; Mira variables with HBB tend to me more 
luminous than their period would suggest, so the stan- 
dard period luminosity relation over estimates the dis- 
tance. Current estimates in the literature range from 
1.9 kpc (Kolotilov et al. 1998) to 2.7 kpc (Munari et al. 
1990). Throughout this work, we will adopt the Munari 
et al. distance. 

The orbital parameters of the system are not con- 
strained, although Munari et al. (1990) claim an orbital 
period of 43 years based on possible dust obscuration 
events in the historic lightcurve. If we assume that the 
white dwarf mass is ~1 M 0 , and that the Mira mass 
is in the range 4-8 M 0 since it is experiencing HBB 
(Mikolajewska 2010), then the 43 year period implies a 
separation of 20-25 AU. 

V407 Cyg has experienced two small amplitude out- 
bursts (peak visual brightness ~13 mag) in the past, in 
1936 and 1998. The 1936 event was originally charac- 
terized as a classical nova event (Hoffmeister 1949), al- 
though this interpretation seems unlikely given the most 
recent event. The 1998 event has been referred to as 
a symbiotic type outburst. The most recent outburst 
was discovered simultaneously by two Japanese groups 
on 2010 March 10 (Hirosawa 2010; Nakano 2010), with 
a maximum magnitude in unfiltered visual CCD im- 
ages of 6.8 mags. The amplitude of this outburst was 
much larger than the events observed in 1936 and 1998. 
Munari et al. (2011) presented optical and near infrared 
photometry obtained over 200 days following the out- 
burst. From those data, the authors obtained values of 
t 2 and f 3 (the time to fade from peak V band brightness 
by 2 and 3 magnitudes) of 5.9 and 24 days, respectively, 
making this outburst of V407 Cyg a fast nova event. 

Munari et al. (2011) and Shore et al. (2011) also pre- 
sented optical spectroscopy from the first few months of 
the 2010 outburst. The earliest spectra are character- 


ized by bright emission lines of the Balmer series of hy- 
drogen, and of He I. The evolution of these line features 
was almost identical to that observed during the 2006 
outburst of RS Oph. After day 48, a system of coronal 
lines appeared typical of classical novae in outburst and 
again similar to the 2006 outburst of RS Oph. The high 
resolution spectra of Shore et al. (2011) revealed asym- 
metric structure in several emission lines, with red wings 
extending to larger velocities than blue wings. The au- 
thors attribute these features to the asymmetric density 
structure of the binary environment encountered by the 
nova ejecta. The equivalent width of some line features 
(e.g. O I 8448) appears to evolve in tandem with the 
X-ray emission, while others (e.g. O I 6300) more closely 
follow the optical and UV evolution. 

The most remarkable aspect of the 2010 outburst of 
V407 Cyg was the detection of the nova as a transient 
7-ray source with the LAT detector onboard the Fermi 
observatory (Abdo et al. 2010). The first significant de- 
tection of the source was on 2010 March 10, the same 
day as the discovery of the optical outburst. The Fermi 
source is coincident with the optical position of V407 Cyg 
to within 0.040°, well within the 95% confidence posi- 
tional error of the LAT. Both the temporal and spatial 
coincidence with the V407 Cyg optical outburst indicate 
that this event is the origin of the transient y-rays. The 
peak y-ray flux was observed on days 3 and 4 after the 
outburst, after which the source began to fade. The last 
significant (>3er) Fermi detection was on day 15 of the 
outburst (2010 March 25). Gamma-ray emission from 
novae has been predicted at energies <1 MeV due to the 
decay of radioactive isotopes of sodium and beryllium 
(Hernanz & Jose 2005), although this has not observed 
to date. The Fermi detection of V407 Cyg has demon- 
strated that in some cases novae are also capable of pro- 
ducing ~GeV photons, adding a new class of transient 
to the y-ray sky. 

In this article, we use X-ray observations obtained with 
Swift and Suzaku to study the interaction of the nova 
ejecta with the wind of the Mira component in V407 
Cyg. In Section 2, we discuss the available data and our 
reduction procedure. Next, we present a brief optical, 
UV and X-ray overview of the outburst in Section 3 and 
a detailed spectral analysis of the Suzaku data in Section 
4. Then, in Section 5, we analyze the Swift XRT spec- 
tral evolution utilizing the results of our Suzaku analysis. 
In Section 6 we present a simple model of the outburst 
and utilize it to explore the gross features of the X-ray 
lightcurve. We compare the X-ray evolution of V407 Cyg 
with that of RS Oph in Section 7. Finally, we discuss our 
results in Section 8, and present a summary and our main 
conclusions in Section 9. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Swift 

Swift began observing V407 Cyg as part of a target 
of opportunity (ToO) program just three days after the 
outburst. The ToO program was then extended a num- 
ber of times as the system continued to evolve. The final 
dataset consists of 34 individual pointings spanning the 
first 120 days of the outburst. In all but four observa- 
tions, a UV image in the UVM2 filter was also obtained 
with the UVOT. The details of each observation, includ- 
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ing exposure time, XRT count rates and UVOT mag- 
nitudes, are presented in Table 1. An analysis of this 
dataset was also presented by Shore et al. (2011). 

We performed the reduction and analysis of the Swift 
data using a processing script customized for the XRT 
and UVOT data (Donato et al., in preparation). The 
script re-processes the data, stored in the HEASARC 
archive, using the standard HEASoft software (version 
6.8) and the latest calibration database (20091130 for 
XRT and 20100129 for UVOT). The reduction of XRT 
data consists of running xrtpipeline, selecting only 
events with grades 0-12 in photon counting mode (PC), 
while the UVOT event mode data are processed following 
the steps reported in the UVOT Software Guide 2.2. 

The XRT data analysis is performed using the script 
xrtgrblc, available in the HEASoft software package. In 
brief, the script selects the best source and background 
extraction regions based on the source intensity (typi- 
cally a circle and an annulus, respectively). In the case 
of V407 Cyg, the source extraction region has a radius 
of up to 55", and the background annulus has inner and 
outer radii of up to 80" and 135", respectively, when 
the light curve is close to its peak count rate. If held 
sources are present in the background region, they are 
excluded using circular masks whose radius depends on 
the held source intensity. Using these regions, source and 
background count rates, spectra and event lists are ex- 
tracted. The net source count rate has to be corrected for 
irregularities in the exposure map and the Point Spread 
Function (PSF). The total correction factor is obtained 
using the HEASoft tool xrtlccorr. 

The script that handles the UVOT analysis 
(uvotgrblc) determines the presence of held sources and 
excludes them from the background region with circular 
regions whose size again depends on intensity. The 
optimal background region is chosen after comparing 3 
annular regions centered on the main source in summed 
images. The region with the lowest background is 
selected. In the case of V407 Cyg, the held of view 
is contaminated by many CCD artifacts (e.g., ghost 
images and/or readout streaks) from a nearby very 
bright star. For this reason, a custom background 
region has been chosen to minimize the effects of these 
artifacts, consisting of two 15" circular regions located 
40" from V407 Cyg, one at 146° and the other at 326° 
position angle. 

The source extraction region is also intensity depen- 
dent: uvotgrblc selects the aperture size between 3" 
and 5" based on the observed count rate and the pres- 
ence of close held sources. V407 Cyg is a bright source 
even in the UVM2 hlter (A ce „ = 2246 A) and only dim 
held sources are present. A 5" region is selected for every 
pointing. The script estimates the best source position 
using the task uvotcentroid and the photometry using 
the task uvotsource. The resulting values are corrected 
for aperture effects but not for interstellar extinction. 
Magnitudes are presented in the AB system. 

2.2. Suzaku 

We observed V407 Cyg with Suzaku (Mitsuda et al. 
2006) as a target of opportunity between 2010 April 
9 21:32 UT and April 10 19:31 UT (sequence number 
905001010). Here we concentrate on the data obtained 
with the X-ray Imaging Spectrometer (XIS) in the 0.3- 


10 keV band. Of the three active units of the XIS, two 
(XIS0 and XIS3) contain a front-illuminated (FI) CCD 
chip, while one (XIS1) contains a back-illuminated chip 
with higher sensitivity to low energy photons. All three 
were operated in the full- window, imaging mode, obtain- 
ing X-ray event data over the full 19' by 19' field of view 
every 8 seconds. 

After standard screening 9 , there are 42,200 s of good 
on-source data with each XIS. We extracted the source 
event using a 210" radius circular extraction region cen- 
tered on the source. The background was taken from an 
annular region also centered on the source with outer 
and inner radii of 410" and 240", respectively. No 
sources were obvious in the background region. We 
then created the response Hies using the xisrmfgen and 
xissimarf gen tools. Since there were no differences be- 
tween the XIS0 and XIS3 data, we have combined them 
in the subsequent analysis. We ignored data below 0.4 
keV for the combined FI spectrum, and below 0.3 keV 
for the XIS1 spectrum, and grouped the channels so that 
each had at least 25 source counts (in order to facilitate 
the use of x 2 statistics in spectral fitting). The XIS1 data 
appear to have a different normalization, so we have al- 
lowed a cross-normalization term in subsequent fitting. 

We inspected the background-subtracted HXD/PIN 
data to ascertain if V407 Cyg exhibited any excess hard 
emission beyond the extrapolation of the components de- 
tected with the XIS. For the non- X-ray background, we 
used the tuned background model, and we used “typ- 
ical,” high-latitude cosmic X-ray background spectrum 
(Boldt 1987). No signal is detected above 20 keV, with 
a net count rate of —9.7 ± 2.7 x 10 -3 ctss -1 in the 20- 
80 keV range, or —9.6 ± 2.4 x 10 -3 ctss -1 in the more 
sensitive 20-50 keV range. In the 15-20 keV range, an 
apparent signal of +5.7±2.1 x 10 -3 ctss -1 was detected, 
although given the systematic uncertainties in the diffuse 
background, we consider this a non-detection. 

3. AN OVERVIEW OF THE OUTBURST EVOLUTION 
3.1. Optical and UV lightcurves 

The optical and ultraviolet flux from V407 Cyg faded 
over the entire observation period covered in this work. 

The optical V band lightcurve (comprised of verified ob- 
servations obtained from the AAVSO archive) and the 
UV data from the Swift UVOT are presented in Fig. 

1. V407 Cyg had an observed brightness of ~8.5 magni- 
tudes in the V band, and ~12.6 magnitudes in the UVM2 
bandpass at the time of the first Swift ToO. The source 
then faded in all subsequent observations. The rate of de- 
cline was not constant, and not the same in the optical as 
in the UV. From day 5 until day 45 the optical UV emis- 
sion faded slowly, with a decrease of only one magnitude 
in the UVM2 filter. This decay was much slower than 
that observed in the optical over the same time period, 
where the V band magnitude decreased from ~8.5 to ~11 
(see lower panel of Fig 1) . After day 45, the UV bright- 
ness began to decline much more rapidly, decreasing by 
one magnitude in just 12 days . This change in decline 
rate was also observed in the optical lightcurve, and cor- 
responds to a change in the slope of the X-ray lightcurve 
at approximately the same time (see Fig 2.). Finally, 

9 see http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/processing/criteria_xi 
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Table 1 

Summary of observations 


Date Days from to a texp (s) 


iaic ■ ' 


Swift 


2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 

2010 - 


03-13 

03-15 

03-19 

03-23 

03-24 

03-25 

03-27 

03- 31 

04- 02 
04-04 
04-06 
04-08 
04-10 
04-12 
04-14 
04-16 
04-17 
04-18 
04-20 
04-22 
04-24 

04- 28 

05- 02 
05-06 
05-18 

05- 26 

06- 03 
06-11 
06-19 

06- 27 

07- 05 
07-13 
07-21 
07-29 


3.8 

5.8 
9.3 

13.2 

14.7 
16.0 

17.2 

21.2 

23.8 

25.3 

27.7 

29.1 

31.2 

33.1 

35.1 

37.5 

38.1 

39.4 

41.0 

43.8 

45.6 

49.3 

53.3 

57.0 

69.5 

77.4 

86.0 
93.0 

101.9 
109.2 
117.7 

125.9 
133.1 
142.0 


949 

955 
922 

956 
2169 

916 

976 

1237 

3164 

4194 

4217 

3981 

3871 

3002 

2003 

4801 

9828 

4740 

3357 

3602 

2784 

2725 

3266 

2930 

530 

850 

1073 

961 

397 

1236 

1230 

348 

1329 

1140 


0.011 

0.011 

0.023 

0.009 

0.016 

0.023 

0.048 

0.158 

0.210 

0.209 

0.223 

0.244 

0.263 

0.244 

0.235 

0.225 

0.219 

0.218 

0.211 

0.173 

0.149 

0.128 

0.105 

0.089 

0.082 

0.078 

0.091 

0.070 

0.040 

0.046 

0.056 

0.066 

0.033 

0.033 


± 0.004 


0.003 

0.005 

0.003 

0.003 

0.005 

0.007 

0.011 

0.008 

0.007 

0.007 

0.008 

0.008 

0.009 

0.011 

0.007 

0.005 

0.007 

0.008 

0.007 

0.007 

0.007 

0.006 

0.005 

0.014 

0.011 

0.010 

0.009 

0.011 

0.012 

0.008 

0.016 

0.006 

0.006 


<0.012 

<0.008 

0.008 ± 0.003 
<0.014 


0.004 

0.007 

0.017 

0.046 

0.079 

0.089 

0.108 

0.125 

0.161 

0.146 

0.136 

0.150 

0.157 

0.158 

0.148 

0.126 

0.108 

0.095 

0.084 

0.067 

0.056 

0.065 

0.074 

0.066 

0.034 

0.032 

0.045 

0.050 

0.027 

0.029 


0.001 

0.003 

0.004 

0.006 

0.005 

0.005 

0.005 

0.006 

0.006 

0.007 

0.008 

0.006 

0.004 

0.006 

0.007 

0.006 

0.006 

0.006 

0.005 

0.005 

0.012 

0.010 

0.009 

0.009 

0.011 

0.010 

0.007 

0.014 

0.005 

0.006 


0.008 ± 0.003 
0.010 ± 0.003 
0.015 ± 0.004 
<0.009 

0.013 ± 0.002 
0.016 ± 0.004 
0.031 ± 0.006 
0.113 ± 0.010 
0.130 ± 0.006 
0.121 ± 0.005 
0.114 ± 0.005 
0.120 ± 0.005 
0.104 ± 0.005 
0.097 ± 0.006 
0.108 ± 0.007 
0.075 ± 0.004 
0.063 ± 0.003 
0.062 ± 0.004 
0.062 ± 0.004 
0.046 ± 0.004 
0.041 ± 0.004 
0.034 ± 0.004 
0.021 ± 0.003 
0.022 ± 0.003 
0.026 ± 0.008 
0.014 ± 0.005 
0.017 ± 0.004 
<0.011 
<0.021 

0.014 ± 0.007 
0.011 ± 0.003 
<0.030 
0.006 ± 0.002 
<0.008 


12.57 ± 0.03 
12.71 ± 0.03 
12.76 ± 0.03 


12.78 

12.82 

12.90 

12.94 

13.04 

13.03 

13.06 

13.21 

13.21 

13.22 

13.26 

13.30 

13.50 

13.52 

13.86 
14.36 
14.76 
15.41 
15.70 
16.02 

16.15 
16.46 
16.54 
16.72 

16.87 
16.93 

17.15 


0.03 

0.03 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.03 

0.03 

0.03 

0.03 

0.02 

0.03 

0.03 

0.05 

0.03 

0.03 

0.06 

0.04 

0.04 


Suzaku 


2010-04-09 


30 


42200 (FI) 
(BI) 


0.557 ± 0.004 
0.479 ± 0.003 b 


0.300 ± 0.003 
0.227 ± 0.002 b 


0.252 ± 0.003 
0.251 ± 0.002 


Days from 2010 March 10 

^ For Suzaku FI chips (XISO, XIS3), the lower energy bound is 0.4 keV 

around day 60 the rate of decline slowed again, this time 
fading by one magnitude in approximately twenty days. 
Again, this change in the decline rate coincided with a 
similar change in the optical lightcurve, and also with a 
change in X-ray behavior, in this case the cessation of 
X-ray fading and transition to approximately constant 
brightness at day 70. 


3.2. X-ray lightcurve 

The evolution of the X-ray emission over the course 
of the outburst was dramatically different from that ob- 
served in the optical and UV filters. We present the 
0.3-10 keV lightcurve observed with Swift in the upper 
panel of Fig. 2, and the hardness ratio (defined as the 
ratio of the count rates in the 0.3-2 keV and 2-10 keV 
bands) in the lower panel. V407 Cyg was detected from 
the first observation onwards, and the source was initially 
faint (0.011 cts s -1 in the 0.3-10 keV range) and rather 
hard, with twice as many counts above 2 keV as below. 
Over the first two weeks of the outburst there was little 
evolution in either the count rate or the hardness ratio. 

Then, at around day 16, the count rate started to rise 



60 80 100 
Days since outburst 


Figure 1. Upper panel: Swift UVOT UVM2 filter light curve. 
Lower panel: AAVSO V band verified light curve over the same 
time period. 


dramatically, increasing by a factor of 10 from the initial 
values on day 23. The system reached a peak count rate 
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Figure 2. Upper panel: XRT lightcurve in 0.3—10 keV range over 
the course of outburst. Lower panel: Hardness ratio evolution 
with time. The hardness ratio is defined as the ratio of the 2-10 
keV and 0.3-2 keV count rates. 

of 0.26 cts s -1 on day 31 (2010 April 10). As the source 
increased in brightness, the hardness ratio declined, in- 
dicating that the emission was becoming softer. After 
day 31, the emission began to fade, and by day 49 the 
count rate had decreased by a factor of two from its peak 
value, reaching 0.12 cts s -1 . From this point forward, the 
hardness ratio remained approximately flat, with a value 
of ~0.4. The fading continued until day 70, at which 
point the 0.3-10 keV count rate remained approximately 
constant (with some fluctuations) over the subsequent 6 
weeks. 


sorbed only by the ISM (model 1 in Table 2) cannot 
reproduce the data — while the 1-10 keV continuum is 
fitted reasonably well by this model, large residuals are 
seen at lower energies. Adding a second temperature 
component (model 2) does not improve the fit. However, 
the addition of an additional tbabs absorber attenuat- 
ing both plasma components (model 3) results in a dra- 
matically improved fit (x 2 / u = 1-22 for 1087 degrees of 
freedom.) The additional absorber has a column density 
of 7.6 x 10 21 cm -2 , and the two plasma temperatures 
are 0.08 and 3.2 keV. The presence of this additional ab- 
sorber is reasonable given the dense wind of the Mira 
companion. We also tried a blackbody to account for 
the soft flux (model 4). The fit is statistically just as 
good (x 2 /^ = 1-23 for /nu = 1087), and returns similar 
parameter values for the intrinsic absorber (N(H) = 7.2 
x 10 21 cm~ 2 ) and the bremsstrahlung component (kT = 
3.3 keV). The blackbody component has a temperature 
of 57 eV. 

In Figure 3, we show the complete 0.3-10 keV spec- 
trum overplotted with model 4 from Table 2. The resid- 
uals of the model fit are plotted in units of Ax 2 in the 
lower panel of the figure, and clearly demonstrate the 
presence of the emission lines noted above. The most 
significant feature in the model residuals is the He-like 
Fe emission line at 6.7 keV. Since the analysis of the con- 
tinuum showed that the X-ray spectrum is comprised of 
two distinct emission components (dominating the emis- 
sion above and below 1 keV, respectively,) we consider 
each separately as we move forward in the development 
of our spectral model. 


4. A DETAILED SNAPSHOT OBTAINED AT PEAK X-RAY 
BRIGHTNESS 

The Suzaku spectrum, obtained at the peak X-ray 

brightness of the outburst, offers a highly detailed snap- 

shot of the emission of the system. The spectrum is dom- 

inated by continuum emission, but some lines are also 
present, most notably those of He-like Si, S and Fe at 1.9, 
2.5 and 6.7 keV, respectively. There also a distinct soft 
component, peaking at 0.5 keV. Our goal is to develop 
a physically motivated model which reproduces the ob- 
served spectrum. Since the emission is clearly complex, 
we approached this task in stages by focusing on specific 
features of the spectrum and then adding components as 
needed. 

To begin, we ignored all obvious emission features and 
modeled only the continuum. We included a constant 
to account for the different normalization of the FI and 
BI chips, and then tied all model components between 
the front and back illuminated CCDs. In all models, 
we include a tbabs 4 * * * * * 10 neutral absorber using the abun- 
dances of Wilms et al. (2000) to represent the interstel- 
lar absorption, and fix this value at 2.65 x 10 21 cm -2 . 
This value was obtained using the E(B-V) value of 0.5 ± 
0.05 (Shore et al. 2011) derived from various diffuse in- 
terstellar band (DIB) features and the E(B-V) to N(H) 
conversion relation of Predehl & Schmitt (1995). The 
parameters of each model are shown in Table 2, in the 
’’Continuum only” section. 

A single temperature bremsstrahlung component ab- 



Figure 3. Suzaku front illuminated XIS0,3 spectrum (black) and 
back-illuminated XIS1 spectrum (red) obtained on 2010 April 9, 
or day 30 of the outburst, shown with an absorbed blackbody + 
bremsstrahlung model (model 4 in Table 2). This fit was deter- 
mined using the continuum only (see Section 4.1). The fit residuals 
are shown in the lower panel, clearly demonstrating the existence 
of emission lines in the data. 


4.1. The hard (1-10 keV) emission 

For the harder component, we consider only the 1-10 
keV region of the spectrum. In order to model the contin- 
uum and the emission lines at the same time, we use the 
apec thermal plasma models in Xspec. These are models 
of collisional ionization equilibrium (CIE) plasmas devel- 
10 see http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/oped by Smith et al. (2001), and include contributions to 
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the emission from free-free, recombination and line pro- 
cesses. The components of these models are presented in 
Table 2 in the “Hard X-rays only” section. 

We first fit an absorbed single temperature apec model 
to the 1 -10 keV data, using the temperature and absorb- 
ing column found from the continuum fitting as start- 
ing parameter values, and allowing the abundance of the 
plasma to vary freely. We again included a fixed value 
absorber with N(H) = 2.65 x 10 21 cm' 2 in order to rep- 
resent the interstellar absorption contribution. The re- 
sulting model (model 5 in Table 2) has y 2 / 1 ' = 1.44/1449, 
and similar parameter values to the bremsstrahlung fit 
to the continuum (kT = 2.8 ± 0.1 keV; N(H)i„t = 6.9 
± 0.3 x 10 21 cm -2 , where the uncertainties are the 90% 
confidence intervals). The best fit abundance value is 
0.36 ± 0.04 Z q. Examining the residuals of this model, 
the continuum emission and He- like Fe line at 6.7 keV 
are reproduced very closely. However, the large (~4 a) 
residuals at the He-like Si and S lines at 1.9 and 2.5 keV 
that were seen in the simple model fit in Figure 3 remain. 

The existence of these strong He-like lines is at odds 
with the temperature derived for the plasma. In a 2.8 
keV (3.4 x 10 7 K) plasma less than 20% of S and no 
Si at all would be helium- like. In fact, we expect both 
elements to be almost completely ionized at this tem- 
perature if the gas is in collisional ionization equilibrium 
(see e.g. Bryans et al. 2009). We cannot account for 
the presence of these lines by adding additional temper- 
ature components, or multi-temperature plasma models 
such as cemekl. Allowing the abundances to vary does 
not improve the fit either. This is because the plasma 
temperature is determined primarily by the fit to the 
continuum emission, and appears to be well constrained. 
At this temperature, we will never observe He-like lines 
of Si or S no matter how abundant they are, since these 
elements are almost completely ionized. 

Instead, the strength of these lines may indicate that 
some fraction of the emitting gas has only recently been 
shocked, and has not yet had enough time to come into 
CIE. If insufficient time has passed for ionization equilib- 
rium to be achieved by particle collisions, then a larger 
fraction of ions of these species will be in lower ion- 
ization stages, increasing the observed strengths of the 
He-like triplet lines. In order to test this, we added a 
non-equilibrium plasma component to the model, imple- 
mented in Xspec as nei. This model is characterized 
by the temperature, abundance and emission measure 
of the plasma as in the apec models, and also by the 
ionization age r, defined as the product of the plasma 
density n and the time since the passage of the shock 
front, t. This parameter in essence describes how much 
the plasma has progressed towards coming into collisional 
ionization equilibrium since it was shocked. Denser plas- 
mas achieve equilibrium faster since more collisions take 
place. We tied the abundance of the nei component to 
that of the apec plasma, since it contributes the Fe line 
at 6.7 keV, and this provides a strong constraint on the 
abundance of the emitting material. 

The addition of the nei component improves the 
= 1.08/1446), and the new model (model 6 in 
Table 2) accounts for the presence of the He-like lines. 
The intrinsic absorber has a higher value in this model, 
with N(H) = (11.1 ± 0.5) x 10 21 cm -2 . The tempera- 


ture of the nei plasma is 3.9+5r4 keV, with r = (2.4+gJ) 
x 10 10 . The temperature of the apec component is 2.5 
± 0.1 keV, and the best fit abundance is 0.57 ± 0.05 
Z 0 . Comparing the normalizations of the two compo- 
nents, we find that ^80% of the flux originates in the 
apec component, i.e. in plasma that is already in ion- 
ization equilibrium. The remaining 20% is from plasma 
which is still coming into equilibrium, as indicated by 
the ionization age of the nei component. As a check, we 
also tried fitting two nei components. We found that the 
ionization age of one of the nei components tended to- 
wards large values, suggesting a plasma that was already 
in equilibrium. The other nei component had very sim- 
ilar values to the apec + nei fit. The presence of a CIE 
plasma contributing most of the flux appears to be re- 
quired to account for the bright He-like Fe line observed 
at 6.7 keV — at the temperatures found by the fit, Fe will 
not be ionized until r ~ 10 12 cm -3 s (Smith & Hughes 
2010 ). 

The normalizations of the model components are re- 
lated to their emission measures, defined as 

EM = [ n e rn dV ~ n 2 V. (1) 

Jv 

Assuming a distance to V407 Cyg of 2.7 kpc, the emis- 
sion measures of the two model components are 9.9 x 
10 56 cm -3 (the apec plasma) and 2.5 x 10 56 cm -3 (the 
nei plasma). 

4.2. Fitting the entire spectrum 

With a satisfactory model for the hard emission, we 
now expand our analysis to fit the entire 0.3-10 keV spec- 
trum. We learned from the continuum-only fits that the 
soft component could be modeled equally well by ther- 
mal bremsstrahlung emission with kT = 0.08 keV or a 
blackbody with temperature 57 eV. In either case, the 
soft component is subject to the same high line of sight 
absorption as the 1-10 keV emission. We will explore 
both of these possibilities for the soft emission by adding 
to the model developed in Section 4.1. Summaries of the 
resulting models are presented in Table 2 in the “Com- 
plete Spectrum” section. 

We find a reasonable fit to the data (x 2 /^ = 
1.11/1657) fitting the softest emission with a thermal 
bremsstrahlung model (model 7 in Table 2.) The pa- 
rameters of the apec and nei components remain largely 
unchanged, although the model fit implies a higher 
value of the metallicity (0. 8 Z 0 ) than before. The 
bremsstrahlung component has a temperature of 0.05 
keV, slightly lower than that found for the continuum 
fit, and a normalization of 1.4 x 10 , corresponding to 
an emission measure of ~4 x 10 64 cm -3 . The most likely 
origin for such low temperature plasma in this system is 
the reverse shocked ejecta. However, the high emission 
measure seems to rule out this possibility, as it implies a 
very high ejecta mass, M e j. 

To estimate M e j, we assume that the ejecta are com- 
prised only of hydrogen, and that they have been com- 
pletely ionized by the reverse shock, which has already 
reached the white dwarf surface. We further assume that 
the shocked ejecta occupy a roughly spherical volume of 
radius R = v e jt , where v e j is the initial ejecta velocity 
and t is the time since outburst. Taking an ejecta ex- 
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pansion velocity of 3200 km s _1 for 30 days (based on 
the FWZI of the Balmer lines, see Munari et al. 2011 
and Shore et al. 2011), then the emitting volume is ~2 
x 10 45 cm 3 . Inserting this volume estimate into Equa- 
tion 1, we obtain an estimate of the ejected mass of ~7 
x 10 -3 M 0 . This mass estimate can be lowered some- 
what by accounting for deceleration of the ejecta, which 
would reduce the volume. However, even considering a 
volume 10 times smaller than the one assumed here only 
reduces the mass estimate to 2 x 10 -3 M 0 . Such high 
ejected masses would be unprecedented in a nova out- 
burst, and appears to be at odds with the fast evolution 
of the optical lightcurve. 

To allow for the possibility of more complex soft emis- 
sion, we fit the data using an apec model (model 8). 
Although the resulting emission measure is smaller, it is 
still extremely large (a few 10 63 ), implying an ejected 
mass of at least a few 10~ 4 M 0 . Furthermore, the fit 
is not as good (x 2 / u = 1-20/1657), and clear residuals 
are seen in the data at the location of spectral lines of 
e.g. O VIII Ly a at ~0.8 keV. We note that Shore et 
al. achieved acceptable fits to the soft component ob- 
served in the summed Swift spectrum using a variable 
abundance apec model with extremely super-solar N, O 
and Ne abundances. We tried using a similar model to 
fit the soft component in the Suzaku data, allowing the 
abundances of N, O and Ne in the low temperature apec 
component to vary. However, we found no combination 
of parameters resulting in y 2 t v < 1-4, and large residu- 
als are still apparent in the data, particularly between 0.7 
and 1 keV. The emission measure required to reproduce 
the observed flux in this model is lower, although it still 
implies a high ejecta mass at odds with the fast optical 
evolution of the outburst. This contradiction, in combi- 
nation with the extremely high abundances required and 
the poorer fit to the data, lead us to conclude that a 
thermal, line emitting plasma is an unlikely origin of the 
soft emission. 

Finally, we considered a blackbody origin for the soft 
emission (model 9 in Table 2). Including such a source 
is reasonable if we are seeing the supersoft emission from 
the still burning white dwarf. Although the intrinsic ab- 
sorption of the red giant wind is high, it is feasible that 
the supersoft emission is observable due to its high lu- 
minosity (typically > 10 3 ' erg s _1 ). The fit found for 
this model is good, with y 2 / ix = 1.11/1661. Again, the 
parameters for the hard components remain largely un- 
changed from the 1-10 keV only fits. The temperature 
of the blackbody is 38 eV, corresponding to ~440,000 
K. The normalization corresponds to a bolometric lu- 
minosity of 8.4 x 10 39 erg s _1 . While this is unreal- 
istically large, it is now well known from X-ray studies 
of novae during the supersoft phase that blackbody fits 
tend to under-predict the temperature and overestimate 
the luminosity of these shell burning white dwarfs (see 
e.g Osborne et al. 2011 (RS Oph) and Page et al. 2010 
(V2491 Cyg)). The data and this best fit model are pre- 
sented in Figure 4. The 0.3-10 keV flux of the hard com- 
ponent is 9.1 x 10 -12 erg s _1 cm -1 , and the luminosity 
(corrected for absorption) in the same energy range is 2.8 
x 10 34 erg s _1 . 

We made one further refinement to this model. Since 
the intrinsic absorption likely arises in the wind, we tied 
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Figure 4. XIS 0.3-10 keV spectrum with best fit model 9 from 
Table 2. This model has a blackbody for the soft component. The 
lower panel shows the fit residuals in units of sigma. The emission 
line residuals seen in Figure 3 are now gone. 

the abundances of this absorber to that of the apec com- 
ponent, which itself originates in the wind shock heated 
by the ejecta. The fit is statistically the same as the so- 
lar abundance absorber (x 2 /v = 1.11 for v = 1661), and 
all model parameters agree within the uncertainties, sug- 
gesting that the fit is insensitive to the abundance of the 
material, which is not surprising given the large column 
density attenuating the soft flux. 

As an additional check, we extrapolated our best fit 
model into the HXD energies. Using PIMMS v4.3, we 
find that the model predicts a signal of only 2.9 x 10~ 4 
cts s _1 in the 15-20 keV range. This is consistent with 
the non-detection in the HXD detector. 

5. EVOLUTION OF THE TWO X-RAY EMITTING 
COMPONENTS OVER THE OUTBURST 

The model we have developed for the Suzaku data gives 
a detailed snapshot of the interaction of the nova ejecta 
and the Mira wind. We now utilize that model to in- 
terpret the Swift lightcurve, focusing on the evolution 
of the two distinct spectral components observed in the 
Suzaku spectrum. As we discussed in Section 3.1, most of 
the Swift observations were too short to collect enough 
counts for spectroscopy. However, we can increase the 
statistics by summing neighboring observations together 
in order to investigate the origins of the changes in the 
X-ray count rate and hardness ratio we observed in the 
lightcurve. 

During the first two weeks of the outburst the count 
rate remained low with a relatively constant hardness ra- 
tio, ranging between 2 and 3. Since we observed very 
little evolution in the lightcurve over this period, we 
summed the first seven observations (up to March 27) 
together. From this date and forward, we summed three 
consecutive observations at a time to derive spectral in- 
formation (e.g., from March 27 to April 2, April 4 to 
April 8, and so on). If the combined spectrum had less 
than ~200 counts, we fit our model to the unbinned data 
and utilized the Cash statistic to assess confidence in- 
tervals for our parameters of interest. If the resulting 
spectrum had more counts, we binned the data to give 
a minimum of 20 counts per bin, and adopted the x 2 
statistic to assess best fit parameter values and associ- 
ated confidence intervals. Hereafter, when we mention 
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the number of days after the 7-ray trigger, we refer to 
the time of the middle observation in each summed group 
of three. 

The nei component in our best fit Suzaku model re- 
quires good statistics in the emission lines in order to pin- 
point the ionization age and the contribution of this com- 
ponent versus the CIE plasma. Since the Swift spectra 
have rather low numbers of counts (even after summing 
consecutive observations), and since the nei component 
only contributes ~20% of the flux at the time of peak 
brightness, we decided to use a single apec model only 
to approximate the hard flux. Therefore, our Swift model 
takes the form tbabs*tbabs* (bb+apec) . We fixed the 
abundance of the apec component at 0.8 Z@, the N(H) 
of the first absorber to 2.65 x 10 21 cm -2 (again to rep- 
resent the ISM contribution), and allowed the intrinsic 
absorption, temperatures and normalizations of the two 
components to vary. Furthermore, we fixed the abun- 
dance of the intrinsic absorber to that of the apec com- 
ponent, since both the emission and absorption arise in 
the Mira wind. 

The bb+apec model provides a reasonable fit to the 
data for all summed spectra after day 17. The resulting 
parameter values for the two emitting components, the 
total flux and contribution from each, and the column 
density of the intrinsic absorber are plotted in Figure 5. 
The model does not, however, result in a good fit to the 
summed spectrum obtained for the first two weeks of the 
observation (see Figure 6). Although the 2-10 keV flux 
in this spectrum is reproduced well, the model underesti- 
mates the counts at E < 1 keV. There is a large residual 
around 0.4 keV, which seems to arise primarily in the 
earliest observations. It appears that a single blackbody 
cannot adequately account for the flux between 0.3 and 1 
keV. Furthermore, the normalization of the best fit black- 
body returned by the model fitting implies a luminosity 
in excess of 10 43 erg s -1 — a non-physical quantity for this 
system. Although this may be an affect of averaging over 
a relatively large period of time (two weeks), there are 
insufficient counts to draw any strong conclusions about 
the poor fit. 

The parameters of the spectral fits show that the 
shocked gas cools significantly from its highest point dur- 
ing the first two weeks, when the recorded temperature 
was ^6 keV, to the beginning of the plateau around day 
60, when the temperature decreased to ~1.5 keV. After 
this date, the temperature remains approximately con- 
stant. The emission measure of this component peaks 
on day 30 at a maximum value of 1.15 x 10 57 cm -3 , 
in good agreement with the Suzaku model value (1.24 
x 10 57 cm -3 .) The emission measure then drops by al- 
most an order of magnitude over the subsequent 30 days, 
coinciding with the drop in temperature. The decline be- 
comes shallower after day 50, decreasing by only a factor 
of 3 by day 140. 

The blackbody component is first detected after day 
20, although the normalization is poorly constrained un- 
til day 30. Between days 20 and 40, there is an apparent 
increase in the temperature, from 0.03 to 0.06 keV, fol- 
lowed by a drop around day 50 to 0.045 keV. Over this 
time period, the normalization varies between 0.02 and 
1.5, although with large uncertainty, implying bolomet- 
ric luminosities in the range 10 36 to 2.7 xlO 37 erg s -1 . 


As a check, we modeled the individual Swift spectrum 
obtained on day 29, just prior to the Suzaku observation. 
We find that the blackbody temperatures agree closely. 
The blackbody normalization is about an order of mag- 
nitude smaller in the Swift best fit model than in the 
Suzaku model, although within the rather large uncer- 
tainties they agree. The differences are likely due to the 
sensitivity of the blackbody flux to the intrinsic absorp- 
tion in the model — small changes in N(H) correspond 
to large changes in blackbody normalization. After day 
50, the blackbody component is no longer required in 
the model fit, and its normalization becomes consistent 
with zero in all subsequent models. The column density 
of the intrinsic absorber shows a clear downwards trend 
with time, from a maximum value of ^3.2 x 10 22 cm -2 
to zero around day 60. In model fits after this date, 
the absorption is consistent with an origin in the ISM 
only. We note that as the intrinsic absorber decreases in 
strength, the soft emission would become easier to ob- 
serve if present at the same flux level, and so the turn off 
of the blackbody component appears to be real. 

In the lower right of Figure 5, we plot the observed 
fluxes in the 0.3-10 keV energy range of the total model 
(in black), the harder apec (in green) and the blackbody 
(in red). These have not been corrected for absorption. 
The majority of the flux at all stages of the evolution 
(and all flux after day 60) originates in the hard compo- 
nent. The black body can be seen to brighten and then 
turn off. This allows us to understand the hardness ra- 
tio evolution shown in Figure 1 — the evolution towards 
a softer spectral state is only marginally influenced by 
the disappearance of the black body component, and is 
primarily due to the drop in temperature of the hard 
component. 

6. A SIMPLE MODEL OF THE NOVA EJECTA— MIRA 
WIND INTERACTION 

In light of the X-ray evolution observed in RS Oph 
during its 2006 outburst, we propose that the hard X- 
ray component observed in V407 Cyg also originates in 
the forward shock being driven into the AGB star’s wind 
by the nova ejecta. In order to understand the key fea- 
tures of this evolution — the initial low count rate, the 
dramatic rise in brightness around day 20 and the subse- 
quent fading — we have developed a simple model of the 
interaction of the nova ejecta with their environment. 
The model tracks the sweeping up of the Mira wind by 
the ejecta as they travel outwards from the white dwarf, 
and estimates the key properties of the resulting X-ray 
emission: the total emission measure, the weighted aver- 
age X-ray temperature, and the total luminosity. These 
can then be compared to the values obtained from the 
data. 

In spherically symmetric models, where the nova out- 
burst happens at the center of the companion wind, the 
X-ray quantities of interest depend only on the distance 
from the white dwarf since the density of the material 
encountered by the ejecta is the same in all directions. 
This allows the use of simple blast wave approximations, 
as was done for RS Oph (e.g. Sokoloski et al. 2006). How- 
ever, if the white dwarf is offset significantly from the ori- 
gin of the giant companion’s wind, these relationships no 
longer hold since the ambient density depends strongly 
on the direction of motion of the ejecta. The ejecta en- 
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Figure 5. Evolution of the absorbed blackbody + apec model fit to the Swift data. The parameters were determined for 3 day averaged 
spectra. Top row: Evolution of apec emission measure (left) and temperature (center), assuming a distance of 2.7 kpc. Evolution of the 
column density of the intrinsic absorber (right). Bottom row: Evolution of the blackbody normalization and temperature (left and center). 
Evolution of the absorbed flux of the total model fit (black crosses), and the contribution from each component, apec in green diamonds 
and bb in red triangles (right). 
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Figure 6. Best fit of the bb+apec model to the summed data from 
the first two weeks of the outburst. The red line is the combined 
model, and the dashed lines indicate the blackbody and plasma 
components. The data have been rebinned for plotting purposes 
to have a minimum of 5 counts per bin. 


counter a medium of increasing density towards the Mira 
companion, compared with one of decreasing density in 
the opposite direction. 

To account for this asymmetry in the environmental 
density, we utilize the geometry presented in Abdo et al. 



Z7 >i 

A0, 


Figure 7. The geometry of the system used in the model. The co- 
ordinate system is centered at the white dwarf (the white circle), 
embedded within the wind of the red giant, which is located a 
distance a away along the z axis. The ejecta travel outwards from 
the white dwarf. We show one example here of a ring of ejecta, 
with extend A Qi. The discontinuity surface between the ejecta and 
the giant wind is located a distance R^j from the white dwarf, and 
Dij from the giant. 


(2010) and illustrated in Figure 7. The circumbinary en- 
vironment is described using a spherical coordinate sys- 
tem centered on the white dwarf, with the red giant lo- 
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cated a distance a away on the polar axis (where a is 
the binary separation). We assume that the nova ejecta 
are instantaneously ejected as a spherically symmetric 
shell, which we divide into a series of i rings in polar 
angle, each centered on and perpendicular to the polar 
axis joining the two stars. In the j th time step of the 
model run, the i th ring is characterized by its radial dis- 
tance from the white dwarf, Ay , and its polar angle 0y 
from the line joining the two stars. We also take Rij 
to specify the shock discontinuity surface for each ring. 
Since all points on a given ring are equidistant from the 
Mira, they are embedded in ambient medium with the 
same density, which we denote py. For each ring, (9y 
will remain constant in time as we assume radial prop- 
agation of the ejecta away from the white dwarf into a 
medium essentially at rest, and so we denote it 9i moving 
forward. Finally, because ptj is the same for all points 
on a given ring, the rings will remain circular as they 
propagate outwards. 

In our model, we consider i = 180 equally spaced rings, 
each with extent A 8i = 1°. The fraction of the total 
ejected mass in each ring, Mi , is simply the fraction of the 
total volume of the ejected shell that each ring occupies 
at the beginning of the simulation (when the shock wave 
is spherical), and is given by 


= 500 Rq. We set vrg = 30 km s 1 based on the equiv- 
alent width of the Na I D lines at 5890 A observed 
both in quiescence (Tatarnikova et al. 2003b) and out- 
burst (Shore et al. 2011). These values result in po = 4.2 
x 10 -16 g cm~ 3 , or number density no = 4.2 x 10 8 cm' 3 
(assuming a hydrogen atmosphere). 

As the model progresses, we keep track of the velocity 
of each ring, Vjj, the mass of red giant wind swept up 
by each ring rriij (note we will use lower case m to refer 
to the swept up material, and upper case M to refer to 
the ejecta), and the volume that the mass swept up by 
each ring occupies, Vij. Both rriij and Vy are initially 
set to zero. In a given time step, we advance each ring a 
distance A ry = Vijt s t ep , where is the velocity of the 
i th ring and t step is the length of one time step, set to 1 
hour. Therefore, the i th ring will sweep up a mass 

A iriij = p^ / / / r 2 sin[Q)drd0d4>, (6) 

J Rij J Oi JO 

during the j th time step, assuming that the wind density 
encountered by each ring is approximately constant over 
the distance Ary. The total mass swept up by the i th 
ring since the start of the outburst, my , is therefore given 
by 


p0 2 

Mi = M e j / sin(9)d9, (2) 

Je i 

where 9\ and 9i are the two angles defining the an- 
gular extent of the i th ring. The ejecta start at the 
white dwarf surface, and are given an initial velocity, 
uo = 3200 km s _1 , based on the FWZI of the hydro- 
gen Balmer lines observed in the earliest phases of the 
outburst (Munari et al. 2011; Shore et al. 2011). There- 
fore, each ring has an initial reservoir of kinetic energy 
KEi = 0.5 Mivl. 

Assuming that the Mira wind can be described by an 
r~ 2 density profile, the ambient density encountered by 
the i th ring in the j th time step, py can be expressed as 


Pij — Po 



(3) 


my = m,y - 1 + Amy. (7) 

As the rings of ejecta travels outwards, they each sweep 
up and shock the ambient medium they encounter. In 
doing so, they transfer some of their kinetic energy to the 
swept-up material, and so the rings decelerate over time. 
If we make the simplifying assumption that no mixing of 
the ejecta and the swept up material occurs, then My 
the mass of a given ring of ejecta, remains constant with 
time. At the start of the simulation, we assume that the 
evolution of the ejecta-environment interaction is energy 
conserving , such that 


(Mi + my) v\j = (Mi + mij-i) v 2 j_ t . (8) 

At each time step, the velocity of the i th ring is recalcu- 
lated using the above relationship. The post-shock tem- 
perature of the swept-up gas, Ty , is related to vy by 


where Rrg is the photospheric radius of the red giant, 
po is the density of the Mira wind at Rrg , and Dy is 
the radial distance of the i th ring from the Mira in the 
j th time step. Dy is given by 

Dy = \Ja 2 + R^ - 2aRijCos(9i) - Rrg, (4) 

recalling that Ay is the radial distance of the i th ring 
from the white dwarf. The density at the base of the 
Mira wind, po, is 


Mrg 

4n R rg v RG ’ 


( 5 ) 


where Mrg and vrg are the mass loss rate and termi- 
nal velocity of the red giant wind. The quantities Mrg 
and Rrg are not well known for V407 Cyg, and so we 
take fiducial values of Mrg = 10 -7 M 0 yr _1 and Rrg 
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(9) 


At the earliest stages of the outburst, the swept-up mass 
is much lower than the ejected mass, and therefore lit- 
tle deceleration of the ejecta occurs. This stage can be 
considered the free expansion phase of the outburst evo- 
lution. Once the swept-up mass equals the ejected mass, 
the deceleration becomes much more pronounced, and at 
this point the outburst is in the Sedov phase of evolution. 
We assume that the energy conserving phase lasts until 
the temperature of the post shock swept up gas drops to 
10 6 K, at which point the swept-up material can cool ef- 
ficiently via radiation. We refer to this stage of evolution 
as the radiative phase, and at this point the ejecta-wind 
interaction is no longer energy conserving, but momen- 
tum conserving , such that 
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{Mi + rriij) Vij = {Mi + rriij- 1 ) v itj - 1 . (10) 

Since gas with temperatures below 10 6 K emits very lit- 
tle X-ray flux, we assume that swept-up material in the 
radiative phase does not contribute to the X-ray emis- 
sion. 

In each time step, we calculate the emission measure 
of the material swept up by each ring using Equation 1. 
This material occupies the region between the disconti- 
nuity surface Rij , and the forward shock front, which we 
assume is at a distance of 1.2 Rij from the white dwarf 
(see e.g. Chevalier 1982). The factor of 1.2 is appropri- 
ate for ejecta expanding into anpoc r~ 2 profile medium. 
The associated volume of this region, which we denote 
Vij , is given by 


/* 1.2 Rij nOi r2ir 

Vij = / / r 2 sin{6) dr d6 d(f). (11) 

J Rij J &i- |_i Jo 

Assuming a constant density across the forward shock 
region, the number density of the swept up material is 
given by 


riij = rriij / nm H Vij , ( 12 ) 

recalling that rriij is the total mass swept up by the i th 
ring since the onset of the outburst. This can be inserted 
into Equation 1 to obtain the emission measure: 

EMij = n 2 j Vij . (13) 

Finally, we can combine this emission measure and the 
shock temperature in Equation 9 to obtain an estimate 
of the X-ray luminosity of each swept up ring using the 
following approximation: 

Lij = 1.4 x 10" 27 T°- 5 n e niZ 2 g B V tj 

~ 1.4 x 1CT 27 T% 5 EMij. (14) 

where we have assumed a completely ionized hydrogen 
plasma and g B is the Gaunt factor, ~1 (see Equation 
5.15b, Rybicki & Lightman 1979). 

The equations laid out thus far calculate the shock 
properties for each individual ring of ejecta. We of course 
have no observational capability to spatially resolve the 
ejecta as they interact with the giant wind. Therefore, we 
must estimate global parameters that can be compared 
with the data by combining the parameters from all rings 
at each time step. The total emission measure £ Aij and 
luminosity Cj are given by a simple summation over all 
rings: 


180 

£M,=J2 EMij , (15) 

i= 1 
180 

= ( 16 ) 

1=1 

The observed temperature, Tj, is an average over all 
rings weighted by their emission measure: 


E i=iEMijTij 

J E-=i EM^ 


(17) 


We track these global parameters over the course of the 
model run, and plot these quantities as a function of time. 
Two examples of these output plots are shown is Figure 8. 
We also produce a visualization of the spatial evolution 
of the ejecta by plotting Rij for all 9i at intervals of 100 
hours. We color code each surface depending on what 
phase (free expansion, Sedov or radiative) the ejecta are 
in. Examples of this are given in Figure 9. Note that 
for the runs shown in Fig. 9, the expansion of the ejecta 
never becomes radiative. 


6.1. Modeling results 

We present the results of our fiducial model, assuming 
an ejected mass of either 5 x 10~ 7 M 0 or 5 x 10 -6 M 0 
(typical of ejecta masses in fast novae, see, e.g. Yaron et 
al. 2005), in Figures 8 and 9. We also explored a small 
parameter space of nova ejecta mass, binary separation 
and wind mass loss rate, since none of these parameters 
are tightly constrained. Independent of the assumed bi- 
nary separation, ejecta mass and red giant mass loss rate, 
we see some gross features in our model results which are 
common to all parameter sets. 

The emission measure and luminosity are initially very 
low. As the ejecta emitted close to 9 = 0° approach 
the surface of the giant companion, these parameters in- 
crease rapidly and dramatically to a peak value several 
orders of magnitude higher than at the start of the sim- 
ulation. Then, once all ejecta are moving outward from 
the central binary, the emission measure and luminos- 
ity begin to drop, although at a rate which is decreasing 
with time. 

The temperature on the other hand starts off at a high 
value set by the initial ejecta velocity. While all the ejecta 
are in the free expansion phase, this temperature remains 
approximately constant. However, the weighted mean 
value begins to drop rapidly around the same time as 
the luminosity increases, since the majority of the emis- 
sion comes from gas very close to the red giant, which 
has slowed the ejecta down significantly due to its high 
density. The temperature then stays approximately con- 
stant for the remainder of the simulation. These features 
can be compared with the evolution of the luminosity and 
temperature observed in the Swift data presented in Fig- 
ures 5. For the models presented in Fig 8, we find close 
agreement in terms of the shape of the lightcurve and the 
observed temperatures. The observed temperature drop 
is steeper in the model than the observe in the data dur- 
ing the first 10-20 days. The most obvious discrepancy 
between the model and the data is found in the emission 
measure and luminosity — the model predicts values that 
are an order of magnitude larger than we see in the data. 
This, and the incorrect rate of temperature change, are 
probably due to some oversimplifications in our model 
assumptions that we will discuss in the next section. 

The choice of input binary separation, ejected mass 
and mass loss rate do not alter the overall shape of our 
curves, but do result in significant changes in the on- 
set of the X-ray rise, the peak emission measures and 
luminosities reached, the post shock temperatures, and 
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Figure 8. Example output of the simple model described in Section 6, showing the evolution of the forward shock emission measure (left), 
emission measure weighted temperature ( middle ) and luminosity (right). These two model runs assumed a binary separation of 25 AU and 
a mass loss rate from the giant of 10 — 7 Mq yr — 1 . In the top row, the nova ejected mass is 5 X 10 7 Mq, and in the bottom it is 5 X 10 -6 

M 0 . 
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Figure 9. Example of the evolution of the size of the nova ejecta over the first ~200 days of the outburst, for an ejected mass of 5 x 
10“ 7 Mq (left) and 5 x 10“ 6 Mq (right). Each line corresponds to an increase of 100 hours since the onset of the outburst. The red circle 
represents the giant, and its radius is to scale. Gold lines correspond to regions which are in free expansion, green lines to regions which 
are in the Sedov (energy conserving) phase. The same binary separation and wind mass loss rate are assumed as in Figure 8. 
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the rate of decay from maximum. We find that for bi- 
nary separations less than 20 AU, the peak X-ray flux 
was reached at much earlier times than observed in the 
data, in agreement with the binary separations implied 
by the presence of HBB in the Mira component (see dis- 
cussion in Section 1.2). Figure 8 demonstrates the effect 
of changing the ejected mass for a single choice of binary 
separation and red giant mass loss rate. Increasing the 
ejecta mass from 5 x 10 -7 to 5 x 10~ 6 M 0 results in the 
time of peak X-ray brightness occurring ~1 week earlier. 
This is because the higher initial kinetic energy of the 
ejecta results in less deceleration by the ambient wind, 
and thus a faster crossing of the binary system. Note 
that the peak emission measures are similar, since this 
is primarily set by the density structure of the red giant 
wind. 

There are significant differences in the plasma temper- 
ature at late times, with the value trending towards ^6 
keV for the higher ejecta mass, rather than 2 keV in the 
1CT' M 0 case. Again, this is because the smaller deceler- 
ation of the ejecta by the wind results in larger ejecta ve- 
locities and hence higher post-shock temperatures. The 
higher temperature therefore results in a slightly higher 
X-ray luminosity. A faster decay from maximum is also 
found in the model with higher ejecta mass, presumably 
because the ejecta are moving more quickly as they pass 
the red giant and expand out into the larger environment. 

Figure 9 shows the evolution of the size of the ejecta 
over time for the two ejected mass cases. Gold lines in- 
dicate that the ejecta are still effectively in free expan- 
sion, i.e. they have not yet swept up their own mass. 
Green lines indicate ejecta that are in the Sedov phase 
of evolution. The effect of the higher mass is obvious - 
the ejecta remain in the free expansion phase for longer, 
since it takes more time to sweep up the ejected mass. 
The ejecta are also larger in volume, since they have a 
higher velocity for a longer period of time. 

6.2. Limitations of the model 

Our simple model neglects some important physics. 
First, we only consider the evolution of the forward 
shock, and do not track in any way the reverse shock be- 
ing driven back through the ejecta. Second, we consider 
only a simple strong shock with no mixing of the ejecta 
and the swept up gas. In reality, important instabilities 
form which lead to mixing across the discontinuity sur- 
face, and regions of high density enhancements which can 
result in rapid cooling of portions of the swept up mate- 
rial. Furthermore, we assume that the gas in the swept 
up shell is thoroughly mixed on a time scale much smaller 
than one time step in our model, so that its density is 
uniform and always given by Equation 12. In reality this 
is not the case, and the swept up shell likely has some 
density profile in the radial direction away from the white 
dwarf. Such a profile will have implications for the value 
of the emission measure, and for how cooling proceeds 
within the shell. In fact, it is likely that the very earliest 
material swept up by rings traveling away from the giant 
cools within a time period shorter than our model run. 
Therefore, assigning this material the shock temperature 
for later material is incorrect, and probably biases us to 
larger values of the emission measure and luminosity. 

Thirdly, we assume that the forward shock is located 
at distance 1.2 R^j from the white dwarf. This value is 


taken from Chevalier (1982), and assumes that the ejecta 
are centered at the origin of an r~ 2 density profile. This 
assumption clearly breaks down for material moving to- 
wards the red giant at early times, and in fact for this 
material the forward shock is likely located closer to the 
discontinuity surface. In these directions, our emission 
measures are probably underestimated as the swept up 
material in fact occupies a smaller volume. 

Finally, we neglect the effects of particle acceleration 
on the temperature of the post-shock gas. The detec- 
tion of V407 Cyg with Fermi during the first few weeks 
is clear evidence that some of the shock energy was im- 
parted into relativistic particles. Tatischeff & Hernanz 
(2007) discuss the importance of this effect on the early 
evolution of RS Oph, and suggest that the diffusive shock 
acceleration of particles was responsible for the observed 
rate of deceleration of the ejecta, which was higher than 
predicted by standard shock wave models. The primary 
piece of observational evidence for this effect are X-ray 
temperatures lower than expected from the observed line 
velocities using Equation 9. Thus, our model may over- 
predict the post-shock temperature at early times. If we 
take the presence of GeV emission to indicate that diffu- 
sive shock acceleration is an important source of cooling, 
then we would expect to see lower X-ray temperatures 
during the first two weeks of the outburst. Unfortunately, 
this is the epoch where the X-ray emission is very faint, 
and we cannot determine from the existing data if this 
was an important source of cooling. 

Clearly, the assumptions we have made have implica- 
tions for the parameter values (emission measure, tem- 
perature etc) implied by our simple model. However, we 
note that regardless of the exact details of the physics, 
most of the of X-ray emission in our model originates in 
the material swept up by the ejecta traveling towards the 
red giant. 

7. COMPARISON WITH RS OPH 

To date, RS Oph is the only other nova with a gi- 
ant companion which has been well observed in X-rays 
(see, e.g. Bode et al. 2006; Sokoloski et al. 2006). Here, 
we contrast and compare their X-ray evolution. During 
the 2006 outburst of RS Oph, the Burst Alert Telescope 
(BAT) onboard the Swift satellite obtained a statistically 
significant detection of the system in the 14-25 keV range 
in the first five days. Despite the detection of high energy 
gamma rays with Fermi , no similar BAT detection was 
made during the first few days of the V407 Cyg outburst. 

For both systems, the first Swift XRT observations 
were obtained on day 3 of the outburst, with subsequent 
observations occurring every few days at a similar rate. 
In the case of RS Oph, the first pointing was extremely 
bright, with a 0.3-10 keV XRT count rate of 14.2 ± 0.2 
counts s _1 . Even allowing for a greater distance to V407 
Cyg, this is still several orders of magnitude brighter than 
in V407 Cyg (with 0.011 ± 0.004 counts s _1 on day 3). 
Furthermore, RS Oph reached its maximum 0.3-10 keV 
count rate of 31.5 ± 0.2 around day 5, and then faded 
in all subsequent observations. Given that the similar 
initial ejecta velocities (4000 km s _1 for RS Oph vs 3200 
km s _1 for V407 Cyg) imply shock temperatures that 
differ by only a factor of 1.5, and that the intrinsic col- 
umn densities are also similar, then the large difference 
in X-ray count rate must be due to the very different 
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density of the material in the immediate vicinity of the 
white dwarf. 

The most striking difference between the lightcurves of 
RS Oph and V407 Cyg is the emergence of an extremely 
bright supersoft source in RS Oph (Nelson et al. 2008; 
Ness et al. 2009) , believed to originate in the still hydro- 
gen burning shell on the surface of the white dwarf. In 
RS Oph, this component emerged about 25 days into the 
outburst, and then turned off quite suddenly around 50 
days later. While we do detect a blackbody like com- 
ponent with a similar temperature in V407 Cyg, its ob- 
served count rate is much lower. Fits to the Swift spectra 
presented in Bode et al. (2006) show that the intrinsic 
absorption had dropped to a few 10 21 cm -2 by day 10 
of the 2006 outburst in RS Oph. In stark contrast, the 
intrinsic absorption in our model fits to the V407 Cyg 
Suzaku spectrum on day 30 is ~1.5 x 10 22 cm -2 , almost 
an order of magnitude higher. Therefore, the lack of a 
bright supersoft phase in the outburst evolution of V407 
Cyg is due to the much higher intrinsic absorption and 
soft X-ray attenuation than in RS Oph. 

8. DISCUSSION 

8.1. The role of environmental asymmetry in the X-ray 
evolution of the outburst 

The observed brightening of the lightcurve around day 
30 can be explained by the approach of the ejecta towards 
the red giant photosphere. Even though only a small 
fraction of the ejecta intersects the red giant, the large 
increase in density towards the base of the red giant wind 
results in a dramatic increase in the emission measure of 
the shell, and therefore of the X-ray luminosity. The 
rest of the shell emits at a lower luminosity and higher 
temperature, as deceleration of the ejecta is less efficient 
in these directions. This is reflected in our best fit Suzaku 
model, where we find that ~80% of the 2-10 keV flux 
originates in a 2.5 keV plasma that is already in CIE, and 
the remaining 20% of the flux is due to the hotter NEI 
plasma, with a temperature of ~3.9 keV. We examined 
the output of our model, and found that by dividing the 
emitting shell into two regions contributing 80% and 20% 
of the flux, the emission measure averaged temperatures 
of these two regions were similar to those found in our 
spectral fits. Furthermore, the region emitting 80% of 
the flux was within 20 degrees of the line of sight between 
the two stars, i.e. a small fraction of the total swept up 
shell. 

The presence of both NEI and CEI plasmas is firm ob- 
servational evidence of the direction dependency of the 
density profile encountered by the nova ejecta. Accord- 
ing to Smith & Hughes (2010), most species will be in 
collisional equilibrium with an ionization age r = 10 12 
cm” 3 s. In the vicinity of the red giant, the density is high 
and plasma comes into equilibrium extremely quickly. In 
the opposite direction, the density is lower, and the time 
to come into equilibrium is longer. Also, since the den- 
sity in this direction drops as the ejecta move outwards, 
the equilibrium timescale will decrease with time. The 
presence of NEI plasma over the course of our 40 ks ob- 
servation, and a best fit ionization age of ~2.2 x 10 10 
cm” 3 s, implies a maximum density of 5 x 10 5 cm” 3 for 
the plasma away from the red giant. Assuming the same 
fiducial values for the red giant radius and mass loss rate 


as in Section 6, this density implies a distance from the 
red giant center of ~10 15 cm, or 67 AU. Assuming lit- 
tle deceleration in the direction away from the red giant, 
and an initial ejecta velocity of 3200 km s” 1 , then the 
outermost parts of the ejecta should be about 55 AU 
away from the white dwarf surface. Adding the binary 
separation of at least 20 AU, then ejecta should be ~75 
AU away form the red giant center. Thus the distances 
are in approximate agreement, and it is reasonable that 
we see an NEI component at day 30. 

8.2. A massive white dwarf in Vf07 Cyg 

Several features of the optical, UV and X-ray evolu- 
tion provide evidence that the white dwarf in V407 Cyg 
is very massive. The optical evolution was fast, typical of 
novae on more massive white dwarfs (like the recurrent 
novae). Our model fits to the Swift and Suzaku spec- 
tra indicate blackbody temperatures in the range 40-60 
eV, a range very similar to that observed in RS Oph 
(Osborne et al. 2011) which is known to host a massive 
white dwarf. The optical and X-ray properties of V2491 
Cyg suggest that it may also be a recurrent nova hosting 
a massive white dwarf (Ness et al. 2011; Darnley et al. 
2011). The blackbody temperatures found for the Swift 
spectra of V2491 Cyg were somewhat higher, averaging 
~70 eV (Page et al. 2010). Despite the uncertainties in 
the model fits, and the problems of using blackbodies to 
model novae in outburst, the presence of flux between 
0.5 and 0.8 keV associated with a soft spectral compo- 
nent seem to favor an emitting source with characteristics 
similar to both RS Oph and V2491 Cyg. 

Considering the X-ray and UV lightcurves in tandem, 
we note that the break observed in the optical and UV 
around day 50 coincides with the turning off of the soft 
component in the X-ray spectrum. A similar behavior 
was observed in RS Oph by Hachisu et al. (2007), who 
observed a break in the optical lightcurve (in that case, 
the Stromgren y-band filter) around day 75. At the same 
time, a rapid fading of the supersoft emission began. The 
authors attribute this behavior to the cessation of nu- 
clear burning on the white dwarf surface. If we assume 
that the observed break and corresponding turn off of 
the soft component in V407 Cyg have the same origin, 
then the short turn off time (~50 days) points to a nova 
outburst on a very massive white dwarf. More massive 
white dwarfs have higher surface gravities and therefore 
reach the critical conditions for a surface thermonuclear 
runaway with much smaller shell masses. 

Finally, the presence of Li in the spectrum of the AGB 
component in V407 Cyg is indicative HBB. Mikolajewska 
(2010) points out the HBB only occurs in stars with ini- 
tial masses in the range 4-8 M 0 , implying that the Mira 
variable in V407 Cyg is an intermediate mass star. Since 
the white dwarf is more evolved, its progenitor must have 
been more massive than the Mira companion. Therefore, 
by evolutionary considerations alone, we would expect a 
massive white dwarf in this system. The observed prop- 
erties of the outburst appear to be in agreement with 
this expectation. 

Other symbiotic stars hosting massive white dwarfs 
(e.g. T CrB, RT Cru and CH Cyg) have been detected 
at energies above 10 keV with hard X-ray instruments 
onboard Swift , Integral or Suzaku (see, e.g. Kennea et al. 
2009). In general, their spectra are well described by 
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cooling flow models with maximum temperatures be- 
tween 30 and 70 keV. This emission originates in shocked 
gas in the accretion disk boundary layer, where disk ma- 
terial must shed its excess energy before accreting onto 
the white dwarf surface. The characteristic temperature 
of this emission increases with white dwarf mass. If sim- 
ilar high temperature emission was detected in future X- 
ray observations of the system in quiescence, this would 
be further proof of a high mass system (although we note 
the counter-example of RS Oph, where despite the known 
high mass of the white dwarf, the boundary layer emis- 
sion in quiescence has a temperature of only 6 keV. For 
further details, see the discussion in Nelson et al. (2011)). 

8.3. Why the X-ray rise and j-ray decline are 
temporally coincident 

Abdo et al. (2010) first noted the striking temporal co- 
incidence between the cessation of the 7-ray emission 
observed with Fermi , and the rapid brightening of the 
X-ray flux observed with Swift. Our detailed X-ray anal- 
ysis and simple modeling of the outburst have identified 
the sweeping up of material in the immediate vicinity of 
the red giant as the source of rapid X-ray rise. This is 
the opposite of the suggestion in Abdo et al. that the 
X-rays rise due to “the increasing volume of shocked gas 
in the nova shell expanding in the direction away form 
the red giant”. Knowing that a fraction of the ejecta 
have reached the red giant at the time of the gamma ray 
decline allows us to propose a scenario for the gamma 
ray turn off. 

The GeV gamma-ray emission observed with Fermi is 
produced by the interaction of a population of particles 
accelerated to high energies in the shock front with a 
second group of particles (Abdo et al. 2010). In the lep- 
tonic scenario, accelerated electrons up-scatter infrared 
photons to GeV energies through the inverse Compton 
process. In the hadronic scenario, accelerated protons 
interact with slow protons in the swept up medium, pro- 
ducing neutral pions which subsequently decay and emit 
GeV photons. A detectable signal therefore requires a 
combination of 1) a shock speed that can accelerate par- 
ticles up to energies high enough to produce GeV photons 
through these secondary interactions and 2) a number of 
accelerated particles and subsequent interactions suffi- 
cient to produce a high gamma ray flux. 

In directions away from the giant, the density encoun- 
tered by the ejecta decreases. Therefore, the number of 
accelerated particles will drop with time, even if the shock 
velocity remains high and can still accelerate particles to 
high energies. Furthermore, the density of target parti- 
cles is also dropping off with time for both the leptonic 
and hadronic production scenarios, since the radiation 
field and the Mira wind both have 1/r 2 density profiles. 
In combination, these two effects result in a rapidly de- 
creasing gamma ray flux as the ejecta travel out from the 
white dwarf. 

In directions towards the red giant, the density of the 
wind and the radiation field is increasing, and so the 
number of accelerated particles and targets will grow 
with time. While the shock velocity remains high (i.e. 
while little mass is swept up), the gamma ray flux re- 
mains high. However, as the ejecta approach the red 
giant, they experience a rapidly increasing deceleration 
as more and more red giant wind is swept up. This drop 


in shock velocity reduces the maximum energy that the 
particles can be accelerated to. This will cause the GeV 
flux to decrease rapidly, despite the fact that both the 
number of accelerated particles and targets is increasing. 
As the shock stalls, there is simply insufficient energy 
available to power the production of GeV photons. This 
makes it likely that most of the gamma rays were pro- 
duced in the material swept up by ejecta moving towards 
the red giant. We note that while we have presented 
a plausible explanation for why the coincidence in the 
gamma ray decline and the the X-ray rise based on our 
improved understanding of the X-ray evolution, we can- 
not use this new insight to favor a leptonic production 
mechanism over a hadronic one for the gamma rays. This 
is because the radiation field and the Mira wind have the 
same density profile. 

Given the requirements outlined here for a significant 
gamma ray flux, we expect that only nova outbursts in 
wide binaries, and with giant companions, will ever pro- 
duce gamma ray emission of the type observed in V407 
Cyg. The high density environment provided by the gi- 
ant wind is crucial for supplying the required large num- 
ber of gamma-ray producing interactions. The wide bi- 
nary separation is necessary to allow time for the system 
to be detected - too small a separation and the ejecta will 
decelerate very quickly, making detection with an instru- 
ment like Fermi unlikely. In light of the rarity of sym- 
biotic stars, gamma rays from novae will be extremely 
rare, making the detection of V407 Cyg with Fermi a 
truly extraordinary event. 

9. CONCLUSIONS 

We have presented X-ray observations obtained with 
Swift and Suzaku during the 2010 outburst of the symbi- 
otic Mira V407 Cyg. The X-ray emission is initially faint, 
and then brightens rapidly approximately two weeks af- 
ter the start of the outburst. A detailed spectroscopic 
analysis of the X-ray emission at peak brightness reveals 
the existence of two distinct emitting components. The 
first dominates between 2 and 10 keV, and originates in 
the forward shock being driven into the Mira wind by 
the ejecta. This emission is itself quite complex. We find 
evidence for two different plasmas, one in collisional ion- 
ization equilibrium with kT = 2.3 keV, contributing ~80 
% of the observed flux, and one non-equilibrium ioniza- 
tion plasma with ionization age r = 2.2 x 10 10 cm' 3 
s and kT = 4 keV, contributing the remaining 20% of 
the flux. The presence of an NEI plasma component is 
direct observational evidence of the asymmetry in the 
environment encountered by the nova ejecta. 

The second component peaks at energies below 1 keV, 
and we identify it as supersoft X-ray emission from the 
white dwarf seen through the dense Mira wind. This 
emission appears to turn off around day 50, coinci- 
dent with a change in the slope of the optical and UV 
lightcurves, and therefore most likely indicates the ces- 
sation of nuclear burning on the white dwarf. Such fast 
turn off times are indicative of nova outbursts on massive 
white dwarfs, including the recurrent novae. Therefore, 
V407 Cyg joins the set of symbiotics hosting massive 
white dwarfs. Future observations of a high-temperature 
X-ray spectrum in quiescence will provide an important 
test of this result. 

We have developed a simple model of the ejecta- 
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environment interaction, and have used it to show that 
the rapid rise and subsequent decay in the X-ray light 
curve is due to the ejecta approaching and passing the 
red giant. As a result, the most likely explanation for the 
coincidental gamma-ray decline/X-ray rise is the stalling 
of the shock front as the ejecta approach the red giant 
photosphere. We expect that GeV emission will only be 
detectable in nova events occurring in wide binaries with 
giant companions, and will therefore be extremely rare. 

We thank the Swift and Suzaku projects for the gen- 
erous allocation of target-of-opportunity time that made 
this work possible. We also thank Martin Laming and 
Cara Rakowski at the Naval Research Laboratory for in- 
teresting discussions about shocks. 
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